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Abstract 

Lattice QCD with two flavors of Wilson fermions can exhibit spontaneous breaking of flavor 
and parity, with the resuhing "Aoki phase" characterized by the non-zero expectation value 
{il^^^T^ip) 7^ 0. This phenomenon can be understood using the chiral effective theory appropri- 
ate to the Symanzik effective action. Within this standard analysis, the flavor-singlet pseudoscalar 
expectation value vanishes: {iip'y^ip) = 0. A recent reanalysis has questioned this understand- 
ing, arguing that either the Aoki-phase is unphysical, or that there are additional phases in which 
{iipj^ip) 7^ 0. The reanalysis uses the properties of probability distribution functions for observables 
built of fermion flelds and expansions in terms of the eigenvalues of the hermitian Wilson-Dirac 
operator. Here I show that the standard understanding of the Aoki-phase can, in fact, be con- 
sistent with the approach used in the reanalysis. Furthermore, if one assumes that the standard 
understanding is correct, one can use the methods of the reanalysis to derive lattice generalizations 
of the continuum sum rules of Leutwyler and Smilga. 
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I. INTRODUCTION 



New patterns of spontaneous symmetry breaking can arise in lattice theories away from 
the continuum hmit. This paper concerns the example of lattice QCD (LQCD) with two 
flavors of (possibly improved) Wilson fermions. Aoki proposed long ago that the apparent 
masslessness of the pions at non-zero lattice spacing could be understood if there is a phase 
in which parity and flavor are spontaneously broken p^]. Numerical evidence for such a 
phase (largely in quenched studies) |2] was subsequently supported by theoretical analyses 



based on the linear sigma-model [3| and on applying chiral perturbation theory (xPT) 
to the continuum effective Lagrangian of Symanzik [4]. The latter analysis incorporates 
discretization errors in a systematic and theoretically well-established way, and uses the 
standard methods of continuum xPT. Near the continuum limit, it predicts two possible 
scenarios, depending on the sign of an unknown low-energy coefficient. In one scenario, 
flavor and parity are broken, and there is an Aoki-phase, while in the other (the "flrst- 
order" scenario) there is no spontaneous breaking of lattice symmetries. The present paper 
concerns the former scenario, and thus assumes that there are choices of gauge and fermion 
actions which lead to the appropriate sign of the low-energy coefficient.^ 

This analysis of the Aoki-phase has been recently questioned by Azcoiti, Di Carlo and 
Vaquero (ADV) [8], who study the pattern of spontaneous symmetry breaking (SSB) using 
probability distribution functions (p.d.f.s) of fermionic bilinears Given certain assump- 
tions, they argue that either the pattern of non-vanishing condensates differs from that in 
the usual Aoki-phase (and that hermiticity is violated), or that there are phases additional 
to the Aoki-phase in which there is a differing pattern of condensates. Both possibilities 
are in contradiction with the standard xPT-based analysis. The xPT analysis (recapped 
below) predicts a non-zero value only for the condensate (iip^y^Tsip) , and predicts only a 
single symmetry-breaking phase. 

The main purpose of this paper is to point out a loophole in the assumptions made by 
ADV, one that allows the standard Aoki-phase to be present. There is no need for the 
presence of more exotic phases. Furthermore, if one assumes that the standard analysis is 



^ Recent unquenched simulations with unimproved twisted-mass fermions in fact indicate that the first- 
order scenario apphes for both the Wilson gauge action \^ and the tree-level Symanzik improved gauge 
action Q , although an Aoki-phase is observed with the Wilson gauge action at stronger coupling f?! . 
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correct, one can derive sum-rules for the eigenvalues of the hermitian Wilson-Dirac operator 
on the lattice, analogous to those derived by Leutwyler and Smilga in the continuum flo| . 

It will be useful to have in mind the main players in the following discussion. These are 
the zero-momentum parity and parity-flavor violating bilinears: 

Co = 1^ ^^75^(«) , and C3 = — ^ ^ 27/^75X3^(71) . (1.1) 

-''site n -'''site n 

Here n labels lattice sites, A^site is the number of such sites, and the flelds are bare lattice 
flavor doublets. Note that Co and C3 are dimensionless lattice quantities. 

The remainder of this paper is organized as follows. The next section provides a brief 
summary of the xPT argument of Ref. [3], and gives some additional results for the "e- 
regime" that will be needed later. Section HTTl then gives a brief summary of the p.d.f. method 
and the results relevant for its application to the Aoki-phase. An important point is that 
the results from the p.d.f. analysis are completely consistent with the predictions of xPT, 
so that what is new about the approach of ADV is their use of eigenvalue decompositions. 
In Sec. HV] I summarize the argument of ADV, pointing out its inconsistency with the xPT 
analysis. The core of this paper is Sec. |Vl in which I explain the loophole in the argument 
of ADV, and give two examples of how this loophole might apply. I also present the above- 
mentioned sum-rules. I conclude with a summary and a brief discussion of generalizations 
in Sec. |Vll I include two appendices, the flrst describing the derivation of a result used in 
Sec. IVj the second providing an alternative formulation of the sum-rules using microscopic 
spectral densities. 



II. REVIEW OF xPT ANALYSIS OF AOKI-PHASE 



In this section I recall the essential features of the analysis of Ref. 4] • Vacuum alignment 
for two flavors of Wilson fermions with no twisted mass term is determined at leading order 
in xPT by minimizing the potential^ 



V'^(E) 



^Tr(S + S^) + ^ [Tr(S + E^) 



C2 



(2.1) 



Here S is an SU (2) matrix proportional to the quark condensate, and the coefficients have 
magnitudes Ci ~ mAq^j-, and C2 ~ a^Aq^^j^, with m the physical quark mass and a the lattice 



^ The analysis in this section is in an effective continuum theory, so that fields, masses etc. have their usual 
continuum dimensions. 
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spacing.^ To produce an Aoki-phase, the two terms in the potential must have comparable 
magnitudes, which occurs if m ~ o^Aqqj^. This, together with m <^ Aqcd, is the power- 
counting used in this paper. Note that discretization errors linear in a have been be absorbed 
into the critical mass rric. 

The Aoki-phase scenario occurs when C2 > 0. Recalling that Ci is proportional to m 
(with a positive coefficient of proportionality), one finds that for |m| large enough that 
|ci| > 2c2, the condensate is aligned with m just as in the continuum: Sq = (S) = sign(m)l. 
For smaller values of |m|, however, the potential is minimized by So = exp{i6Qn ■ f) with 
cos^o = Ci/(2c2) and h a unit vector. The SU{2) flavor symmetry S UTiU"^ is then 
broken to the U{1) subgroup with U = exp{i(j)n ■ r). This results in two massless Goldstone 
pions, with the third pion having a mass proportional to a. The direction of the condensate 
can be fixed by adding a source term to the action, and the standard choice is to add the 
canonical twisted-mass term fiiip'y^T^ip, with ip the bare lattice SU{2) doublet. This gives a 
mass to the (now pseudo-) Goldstone pions proportional to ^/JI. One then sends the volume 
V to infinity, followed by /i — 0, so that all pions are massive except in the final limit. This 
results in the condensate being aligned as Sq = cos 6*0 + isin6'or3, with the charged pions 
being massless and the neutral pion having m^o ~ aAq^j-,.^ 

An important point for the subsequent discussion is that taking V oo before /i ^ 
implies that m^±L ^ oo, so that one is in the so-called "p-regime" for the charged pions. 
Were one to take the limits in the other order, i.e. /i — and then V oo, then the 
zero-momentum fluctuations in the charged pion directions would be unsuppressed, and one 
would be in the "e-regime" for these modes [14]. I will distinguish between expectation 
values obtained in these two orders of limits using the subscripts "p" and "e" . While this 
has the advantage of linking the limits to familiar names, one should also keep in mind that 
both p- and e-regimes are defined more generally. In particular, they are defined also for 
finite V — as the regimes in which m^L ^ 1 or ^ 1, respectively, with L the box size. I am 
not making use of the full extent of these regimes, but rather only the single points in each 
regime reached in the limits described above. Note that m^oL ~ aL oo in both regimes, 

^ More explicitly, in terms of the low-energy coefficients in the chiral Lagrangian, ci = 2f^BQm and C2 = 
— IGW'WqO^ , where I use the notation of Ref. ll]. The physical quark mass is m = Z^imQ — mc)/a, 
with uiQ the bare Wilson quark mass and rric the critical mass. 

^ This analysis of course receives corrections from higher order terms in xPT, as described in Refs. [12l.ll3l|. 
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since a is held fixed, so that fluctuations in the neutral pion direction are always suppressed. 
This means that both regimes for the lattice theory differ from the corresponding regimes 
in the continuum. For example, the continuum e-regime has unsuppressed zero-momentum 
modes in all three pion directions. 

Using the leading order relations between quark bilinears and S, it follows that, in the 
Aoki-phase, the parity and flavor-breaking condensates are 

^(C3)p = z(f 5o/2)Tr (rg [Sq-SJ]) =-sm0o2fBo, (2.2) 



a3 



{Co)p = (f 5o/2)zTr (So - Sj) = . (2.3) 



The matching factors and powers of a are needed to convert from lattice to continuum 
normalization. Here / and Bq are the standard low-energy coefficients of xPT, in terms of 
which the continuum condensate with a standard mass term is 

{i'ip)cont = {uu + dd)cont = -sign{m)2f Bq . (2.4) 

On the lattice, only the pseudoscalar condensates can be determined, due to the mixing 
of ipip with the identity operator. Thus it is the results of Eqs. (12. 2 p and (12.31) that are 
pertinent. I emphasize two features of these results: (i) The magnitude of the flavor- non- 
singlet condensate varies as one moves across the Aoki-phase, with the maximum at m = 
Ci = being equal to | ('?/'V')cont| when appropriately normalized; (ii) The vanishing of the 
flavor-singlet pseudoscalar condensate is true for any Sq G SU{2), and is not special to the 
particular vacuum of the Aoki-phase. 

n 

This completes the discussion of the analysis of Ref. [4]. For the following, it will be 
necessary to generalize the results in two ways. The flrst is to determine the zero-momentum 
two-point functions of the bilinears in the p-regime. The results are simple (and can be 
written in a way that avoids Z-factors): 

C!)^ = {Cs)l , (2.5) 

2\ /r^ \2 



^o), = (Co); = 0. (2.6) 

In words, only the disconnected contributions remain when iVsite — > oo. This is because, in 
the p-regime, the connected parts receive contributions from massive intermediate states, and 
thus are localized in space. Although one of the l/Ngite factors is canceled by translations, 
the other remains and causes the contributions to vanish when iVgite —>■ oo. 



The second generalization is to the e-regime. As discussed above, zero-momentum fluc- 
tuations of E in the charged-pion directions are unsuppressed, and one must integrate the 
zero-modes over the Goldstone manifold IJ]. Non-zero-momentum modes can be ignored 
at leading order. The result is that the flavor-parity breaking condensate now vanishes 

^3 



Zp Jt,q(^SU{2)/U{1) ^ 



0. 



(2.7) 



By contrast, the average of the square of the zero- momentum mode does not vanish: 



{CD. = ^ 



{ipB,/2f {Tr (ra 
1) ^ 

sin^o2/'So; 



a 

^Zp ) JT.o(iSU(2)/U{l) 

3 \Zp) 



(2.8) 
(2.9) 



Comparing the final expression to Eq. (12. 5p one sees how averaging over the Goldstone 
manifold reduces the size of the squared condensate by the geometrical factor of 1/3. For 
the flavor- singlet pseudoscalar, however, the corresponding expectation values vanish. 



{Co). = {CD. 



0, 



(2.10) 



because the vacuum manifold is simply the origin. 

It is straightforward to extend these results to higher powers of the C^, with results that 
are quoted and used in Appendix [Bl 



III. PROBABILITY DENSITY FUNCTIONS FOR FERMION BILINEARS 

The ADV analysis of SSB for Wilson fermions uses probability distribution functions for 
fermion bilinears. Here I briefly recall the essential properties of these p.d.f.s. The p.d.f. 
is familiar for scalar field theories, where it is defined by inserting 5{C — X^n 'Pn) in 
the functional integral (considering here a real scalar field). It is related to the constraint 
effective potential, P{C) = exp[—Nsitc^consti.{C)] and is a useful tool for studying symmetry 
breaking. In particular, in the absence of source terms, it is invariant under the symmetries 
of the action. For example, if there is a ^ —0 symmetry, and this is spontaneously broken, 
then P{C) = (1/2) [(5(C - Co) + 5{C + Cq)] when N^-.te -> oo, with Co the magnitude of the 
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expectation value in the presence of an infinitesimal source. If the symmetry is unbroken, 
then P(C) = 6(C). 

The construction of a p.d.f. is generalized to fermion bilinears in Ref. [9(. It is a non- 
trivial result that the resulting p.d.f. can be used in the same way as for scalar field theories, 
and, in particular, as a tool to study SSB. In the present instance the p.d.f.s of interest are 
those for Cq and C3. If one could calculate -P(Co) and PiC^) using the method of Ref. [9| 
then one could deduce whether SSB occurs and the nature of any broken phases. 

ADV make particular use of P{Cq) and P{C?^ in the "Gibbs state", which means here 
that one evaluates them with no twisted-mass source term, but takes the iVsite 00 limit. 
As noted in the previous section, this puts the theory in the e-regime, in which symmetries 
are manifest. Thus the expectation values of Ca, a = 0, 3 vanish, because of parity and 
flavor-parity respectively: 

(a), = Jim [ dCa P{Ca) a = . (3.1) 

Note that this result gives no information about SSB, since it holds irrespective of whether 
the symmetries would spontaneously break were the limits taken in the other order (A^sitc — *■ 
00 followed by the twisted- mass source term vanishing). It is also possible, as suggested by 
ADV, that there are several vacua, each with differing patterns of condensates, which are 
averaged over in the Gibbs state, leading to a vanishing result. 

Higher moments of the p.d.f.s are, however, order parameters for SSB. Following ADV, I 
focus mostly on the second moment, for which one expects 

{CD, = \im [ dCa P{Ca) Cl = /geom X 1(^)^1' • (3.2) 

Here /geom is a non-vanishing geometrical factor which depends on the vacuum manifold. It 
follows from (13.21) that, if there is no SSB, and {Ca)p = 0, then {C^e will vanish. On the 
other hand, if there is SSB and {Ca)p 7^ 0, then (C^)e will be non- vanishing. 

The value of /geom can be determined if the vacuum manifold is known. For a single 
vacuum, /geom = 1, while for the continuous U{1) (complex scalar field), SU{2)/U{1) (Aoki- 
phase) and SU{2) (continuum chiral symmetry breaking) manifolds the factor is /geom = 1/2, 
1/3 and 1/4, respectively. These are simply obtained by averaging the squared projection of 
the field onto a fixed axis over the respective manifolds. More complicated vacuum manifolds 
with disconnected components would lead to other, less simple, values of /geom- 



ADV apply this methodology to the Aoki-phase. If the standard analysis holds, then 
(C'o)p = while (C3)p 7^ 0, and one then finds 

(C3^ = K^^0. (3.3) 
{C'o)e = 0, (3.4) 

The issue in the following is whether these results are correct, and in particular, whether 
they are consistent with expressions in terms of eigenvalues of the hermitian Wilson-Dirac 
operator. If not, the standard understanding of the Aoki-phase must be wrong. This 
holds also for the extension of the above results to higher powers of the Ca, which are 
straightforward to derive, and which are quoted and used in Appendix [Bl 

I close this section by noting that the results of the p.d.f. analysis can also be obtained 
using xPT. This is shown by the consistency of Eqs. (12. 9p and ( I2.10p from the previous section 
with Eqs. (13.31) and (13. 4p . This consistency is, in fact, preordained, because the required 
average over the vacuum manifold is identical in the two approaches. Each approach has 
its strengths and weaknesses: xPT adds the specific prediction for {Cs)p, Eq. (12. 2p . while 
the p.d.f. analysis does not require an expansion in m and a. Still, for small m ~ a^Aq^j-,, 
as considered here, one does not need the p.d.f. methodology to pose the puzzle noted by 
ADV. 



IV. THE ARGUMENT OF ADV 



ADV base their argument on the expressions for the expectation values in terms of the 
(real dimensionless) eigenvalues of the hermitian Wilson-Dirac operator, H\y = 'y^Diy. De- 
noting these eigenvalues by \j, and including a bare lattice twisted mass, /xq, one has 

(CS = {Cf)„.e + {Cf),„„ (4.3) 




(4.6) 

The first result on each hne can be used to obtain the "p-regime" expectation values (i.e. 
Nsite oo and then /xq 0). Taking the hmits in this order means that the spectrum 
of eigenvalues becomes continuous. The second result on each line (i.e. that with fiQ = 0) 
gives, once iVsite — ^ oo, the "e-regime" expectation values. 

Consider first the linear moments, Eqs. (14.11) and (14.21) . Applying a parity transformation 
flips the sign of A and /io, and thus averaging over a configuration and its parity conjugate 
leads to a vanishing (Cq) in both p- and e-regimes.^ This is consistent with the standard 
expectations, Eqs. (12.31) and (13. ip . The prediction of (14. 2 p that {C^)^ = is also consistent 
with Eq. (13. ip . In the p- regime, however, one can obtain a non- vanishing expectation 

(C), = - lim hm -2^) = -2vr(p^(0)), = -27rp(0) . (4.7) 

This is the standard Banks-Casher relation applied to the present context, in which pu{^) 
is the eigenvalue density of Hw per unit (dimensionless) volume on a given configuration 
(which can be defined when A^^sitc —* oo) and p(A) is its average over configurations. For the 
Aoki-phase scenario to hold one must have p(0) ^ 0, and to match with the xPT prediction 
(12. 2p requires 

p(0) = p^I^M!^ . (4.8) 
Zp n 

Thus far the results are conventional and uncontroversial. 

The apparent problems with the Aoki-phase concern the expectation value (Cq). In the 
standard scenario, this is expected to vanish in both the e-regime [Eqs. (I2.10p and (13.40 ] 
and p-regime [Eq. (12. 6p ]. The expressions in terms of eigenvalues — Eqs. (l4.3H4.5p — do not, 
however, vanish under parity averaging, reflecting the fact that Cq is even under parity. The 
issue is whether they vanish when A^sitc oo. 

ADV argue that, assuming the standard properties of C| in the Aoki-phase, (Cq) does 
not vanish, in contradiction to the standard Aoki-phase scenario. They reach this conclusion 
by considering the two possibilities for the behavior of the eigenvalues under parity: 



^ At a 7^ one does not have an index theorem and the corresponding zero-modes related to topology. There 
can be isolated zero-modes on some configurations, as discussed further in the next section, but these will 
introduce an additional suppression of /1q from the fermion determinant, and thus do not contribute when 
/zo — > to either of the linear moments. This holds also for the quadratic moments. 
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1. puW = pc/(— A): the eigenvalue distribution for /iq = becomes even in A on each 
configuration when A^^site oo. This is what happens if averaging over an infinite 
volume effectively includes a parity average, which is what one would expect. In this 
case (CQ)disc vanishes in both regimes since the summand in fl4.4l) is odd in A. ADV 
focus (in their section IV) on the e-regime. Since (C|)disc,e vanishes identically [see 
eq. dUD], it follows that 

Thus if (C|)e is non-vanishing, then so is (Cq)^. This is in manifest contradiction with 
Eqs. (12. 9 p and (12.101) [or equivalently with Eqs. (13.31) and (13.41) ]. and the standard 
picture of the Aoki-phase fails. To obtain consistency with (14.91) ADV postulate the 
presence of an additional phase with non-vanishing (Co)p. 

2. puW 7^ pf/(— A): the eigenvalue distribution is not symmetric on a given configuration, 
even in the A^^site —>■ oo limit. In this part of the argument (section III of their paper) 
ADV keep /iq non-zero (they call it while sending A^sito oo, and thus they are 
working in the p- regime.^ Since there are then no exactly massless particles, it follows, 
as discussed above, that connected two-point functions vanish when A'site — ^ oo. ADV 
conclude that ^ 

{CX = lim (Co^)disc = -4 ( f / dX^] ) , (4.10) 

which, by assumption, is negative definite. This is in contradiction with the vanishing 
result (12. 6p expected in the Aoki-phase.'' 

Having considered both symmetric and asymmetric p[/(A), and finding both in contradic- 
tion with the Aoki-phase, ADV conclude that the standard Aoki-phase scenario must be 
incomplete. 

I emphasize that, if the argument of ADV is correct, then one must conclude that the 
standard xPT-based analysis is incorrect. For small enough m and a^, the xPT analysis 



^ In general, puW depends on fiQ. Here, however, fiQ is infinitesimal, so the flip in sign of /io under parity 
does not affect puW- 

^ ADV expand further on the implications of (Cq)^ < 0, which violates hermiticity in the continuum limit. 
For my purposes it is, however, sufScient to show that the result (j4.10[) is in contradiction with the 
Aoki-phase expectations. 
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unambiguously predicts only one possible pattern of symmetry breaking, in which {Cl)p = 
(Co)e = 0. This is simply in contradiction with both of the possibilities enumerated above, 
and in particular with Eqs. (14.91) and (14.101) . 

V. CONSISTENCY OF THE AOKI-PHASE WITH EIGENVALUE SUMS 

The conclusion of ADV's argument is surprising. One would not expect that the Aoki- 
phase scenario, based as it is on a straightforward application of the methods of effective 
field theory, could be invalidated by simple properties of eigenvalues of H\y, especially since 
these properties do not appear to be in conflict with the assumptions of the xPT analysis. 

In fact, there is a loophole in the argument of ADV. This arises because the summand in 
Eq. (14. 4p is infrared divergent in the e-regime. This allows an asymmetry in the spectrum 
which is sub leading as A^^site — > oo, and which does not contribute in the infrared regulated 
p-regime, to nevertheless contribute to e-regime expressions. In this way it is possible for 
(Co)disc,p = while (CQ)disc,e 7^ 0. Disconnected and connected contributions to {Cq)^, which 
have opposite signs, can then cancel, leading to the desired result {Cq)^ = 0. 

I stress that showing this possibility exists simply demonstrates that one cannot make 
a judgment about the Aoki-phase scenario using the basic properties of the eigenvalues 
of Hw- The Aoki-phase is an allowed option. To determine whether it actually occurs, 
however, requires further input, such as that provided by applying xPT to the Symanzik 
effective Lagrangian. Using this input, which predicts the Aoki-phase (assuming C2 > 0), 
one can turn the ADV argument around and deduce sum-rules that the eigenvalues must 
satisfy. These are analogues of the Leutwyler-Smilga sum-rules for continuum QCD [lo| . 

That the expectation values (CQ)disc,p and (CQ)disc,e can differ is an example of the non- 
commutativity of the fiQ and Agito —>■ oo limits. This non-commutativity is familiar 
from the properties of (C3) noted above [see Eq. (14. 7p and preceding discussion]. Another 
example is the behavior of (CQ)conn = (C'Dconn, and it will be useful to describe this as a 
warm-up exercise before explaining the loophole. 

In the p-regime one expects (Co)conn to vanish because only massive intermediate states 
contribute. To see how this works with eigenvalues, one writes 




(5.1) 
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Here I have used the result that averaging over configurations allows one to define a con- 
tinuous density p(A) prior to taking A^site — > oo, and I am assuming that 1/A^site and /io 
are small enough that p(A) does not depend on them. The key point is that the integral 
on the right-hand-side is finite, so that (CQ)conn,p when NsHe oo due to the overall 
factor of 1/A^sito. The only possible sources of divergence in the integral are the infrared 
and ultraviolet regions. Only the first two terms in the Taylor expansion of p about A = 
can lead to infrared divergences, but their contributions vanish for any non-zero /io- The 
ultraviolet divergence (arising because p oc A'^ for A ^ oAqcd) is regulated on the lattice 
because there is cut-off on eigenvalues, Amax ~ 1.^ 

By contrast, in the e-regime, (CQ)conn has the infrared-divergent summand l/(A^siteAj)^. 
Since p(0) 7^ 0, the low eigenvalues are approximately uniformly distributed with spacing 
AA ~ l/[A^site p(0)] and with A min ~ l/[A'sitc p(0)]. They thus give a non- vanishing con- 
tribution to (Co)conn whcu A^^sitc — > c>o [ouc that cannot be represented as an integral over 
p(A)]. This is qualitatively consistent with the expectation (13. 3p from the p.d.f. analysis. 
To agree quantitatively with the xPT result (12.91) requires that the following sum-rule hold: 



This constrains the distribution of the small eigenvalues of H^. Note that it must hold 
separately for each value of m throughout the Aoki-phase. 

With this warm-up completed I now return to main quantity of interest, (CQ)disc- For the 
Aoki-phase to be consistent this quantity must vanish in the p-regime and cancel (CQ)conn 
in the e-regime. Vanishing in the p-regime requires 



Since the integrand is finite in the infrared for any non- zero po? this relation is satisfied if 
the density is symmetric: 



I stress that the condition (15. 4p concerns only the spectrum in the A^site — 00 limit — indeed, 
it is only in this limit that pu (and thus 6pu) becomes well-defined. 

® That the ultraviolet divergence is subleading in l/7Vsito is as in the continuum, as has been discussed in 
Ref. Q. This result applies also for the other eigenvalue sums considered below. 




(5.2) 




(5.3) 



Spui\)=pu{X)-pu{-\) = 0. 



(5.4) 
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The question then is how (Co)conn,e can be non-zero and cancel with (Co)disc,e, i-e. how 
the sequence of equahties in Eq. (14.91) of ADV's first possibihty can fail. From Eqs. (14.41) . 
(14. 5 p and (15. 2p the requirement is that 

2. 




TT 

K^t/ \ V A, / 3 



2 



lini.^dEf 1 ) = -p(0)^ (5.5) 



This relation must hold throughout the Aoki-phase (and, in fact, outside this phase too, 
where p(0) vanishes). The issue is whether the left-hand-side can be non- vanishing given 
the symmetry property (15.40 . This is possible if there is an asymmetry which, while van- 
ishing when A^^sitc ^ oo, is enhanced by the infrared divergence in the summand so that 
J2j l/(A'sitoAj) does not become / dX pu{X)/ X when A^site —* oo. If it did have this limit then 
the resulting integral would vanish given the symmetry of pu{^) and one would be back to 
the first inconsistency noted by ADV.^ 

To show that it is possible for the sum-rule (15.50 to be satisfied I need to recall the 
properties of the spectrum of Hw- There is a wealth of literature on this topic, and I will 
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use particularly the results and insights from Refs. [15|, ll6|, ll7|, Q, Q, Q- As noted above, 
exact zero- modes of H^y are suppressed by its determinant and are not relevant. What is 
important is that the spectrum is known to be asymmetric on almost all configurations. This 
asymmetry is a remnant of the exact zero-modes of the Dirac operator which are present 
in the continuum limit on topologically non-trivial configurations. On the lattice these 
would-be exact zero modes end up as near-zero modes of Hw and lead to an asymmetry, as 
described in more detail in Appendix |X1 I expect that the typical magnitude of the resulting 
asymmetry (defined on a given configuration as the difference between the number of modes 
with Xj > and Xj < 0) depends on a and the lattice volume V = a^A^site as 

|nasym| ~ ClVV Aqqj) . (5.6) 

There are several arguments which support this parametric dependence. The most simple 
is to note that, in the continuum, the typical number of zero modes scales as l^zerol ~ 



mVpD^O), with Pd(0) ~ ^qcd eigenvalue density (per unit volume) of the Dirac 
operator [lo|. This result is for the p-regime, which is appropriate since the neutral pion 
is in its p-regime on the lattice. Now, in the Aoki-phase, symmetry breaking by mass 



The limit in Eq. (j4.5p gives the principal part, which vanishes if pu is symmetric. 
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terms competes with that from O(a^) discretization effects, the latter being dominant in the 
center of the phase. Thus it is plausible that one can use the continuum formula with the 
replacements m a^Aqcj-, and rizem ^asym, leading to Eq. (15.61) . Further arguments in 
support of this relation are given in Appendix |X1 

Another result needed below is that would-be exact zero-modes are expected to have 
eigenvalues shifted to A ~ (oAqcd)'' by discretization effects. Here q is an unknown power 
that I argue below may be g = 3. 

The result (15. 6p for the spectral asymmetry is consistent with the p-regime requirement 
(15. 4p . because one must divide by A'site to obtain the spectral density: 

[dXSpu{X)=\im 1^ = 0. (5.7) 

To satisfy the e-regime requirement (15. 5p requires consideration of how the asymmetry de- 
pends on A. For illustration assume that there are more positive than negative eigenvalues 
on a particular configuration. Due to eigenvalue repulsion, one expects that the extra eigen- 
values will impact the spectrum over the region < A < (aAqco)''- To satisfy (15.51) this 
impact must be appropriately peaked at small A. I give two examples of how this could 
work. 

• The first is simple: I assume that the spectral asymmetry manifests itself by the 
presence of 0(1) extra small (positive) eigenvalues with Xj ~ l/[Ai'siteP(0)], while the 
remaining eigenvalues giving rise to the asymmetry are distributed in such a way as 
to give a contribution to J^j^/i^sitcXj) which vanishes when Ngite — ^ oo. As shown 
by the second example, this requires the asymmetry to be less peaked than a l/\/X 
singularity. 

With these assumptions, I]^ l/(A^sitoAj) ~ p(0), and so the square of this sum gives 
(C'Q)disc,e ~ p(0)^- This has the correct magnitude to allow the sum-rule (15. 5p to be 
satisfied. 

Note that in this example I am not making particular use of the dependence (15. 6p 
of the asymmetry on a and V (nor of the arguments given in Appendix [X]) . It is, 
however, somewhat artificial to assume an 0(1) delta-function-like contribution to the 
asymmetry. 

• In the second example, the asymmetry is spread continuously over the range < A < 

14 



(aAqcD)'- For finite NsHe, the discrete sum over eigenvalues can be approximated by 
an integral aside from end effects which can be accounted for choosing the limits of 
integration appropriately. In this sense, one can consider 6pu{\) for finite NgUe- The 
behavior I assume is 

dpu{X) r- , (5.8) 

in which the key feature is the 1 / a/A singularity. The other factors are chosen so that 
one obtains the desired spectral asymmetry: 



^asym 



Nsite / dX 5puiX) ~ (aAqcD) visits ~ av^A^co • (5-9) 

t/ 



The sum in f l5.5p can then be approximated by an integral with the lower limit regu- 
lated with a quantity of order Amm ~ l/[^site p(0)]: 

E - ^ (aAQCD)-^/V(0)^/^ - (aAQCD)^-^)/^ • (5.10) 

j A/siteAj Jl/[N,itc p(0)] A 

Here I have used p(0) ~ (oAqcd)^ from Eq. (14. 8p . The key point is that the result 
(15.101) has a non-zero (and non-infinite) limit as Nsite — > oo. Inserting (15.101) into 
Eq. (15. 5p then leads to a non-zero result for (CQ)disc,e- 

In addition, recalling that p(0) oc (aAgco)^, one sees that for the a dependence on 
both sides of Eq. (15. 5p to match requires g = 3. 

While this example is perhaps more realistic than the first, I stress that it only works 
if 6p diverges as 1/ \/X and not for other powers. 

I do not know if either of these examples represents the actual behavior. Presumably 
it should be possible to determine more detailed information on the distributions of low 
eigenvalues, as has been done in the continuum limit using the methods of random matrix 
theory. A small step in this direction has been taken in Ref. 

The description of the loophole given above is somewhat awkward and unsystematic. This 
shortcoming can be addressed in part by formulating the required consistency conditions in 
terms of the microscopic spectral density and correlations. This is done in Appendix [Bl 
Also included in this appendix is some discussion of how the consistency conditions extend 
to higher orders (corresponding to sum-rules with higher overall powers of A~^). 



2l| 
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VI. CONCLUSIONS 



The major aim of this paper has been to show the Aoki-phase scenario is not ruled out by 
the arguments of ADV. The examples presented in the previous section demonstrate this — 
the set of possibilities considered by ADV is incomplete. The "survival" of the Aoki-phase is 
consistent with the intuition that one cannot rule out the results of xPT using only general 
properties of eigenvalues of H\y. 

If one accepts the standard analysis of the Aoki-phase, then one finds non-trivial con- 
ditions that must be obeyed by the eigenvalues of the hermitian Wilson-Dirac operators. 
These are the sum- rules fl5.2l) and fl5.5l) . which can also be formulated as constraints on 
integrals of microscopic spectral correlators [Eqs. fIBSp and flB9p respectively]. There are, 
in fact, an infinite set of these sum-rules, involving products of any even number of inverse- 
eigenvalues. While at first sight it may seem daunting that the eigenvalues of H]y must be 
distributed so as to satisfy all these sum-rules, I show in Appendix [B] how each sum-rule con- 
strains an essentially independent eigenvalue correlation function, making it more plausible 
that they can all be satisfied. The need to satisfy an infinite set of sum-rules is not special 
to the lattice theory. Indeed, for actions with chiral symmetry (as in the continuum analysis 
of Leutwyler and Smilga), the eigenvalues of the Dirac operator in each topological charge 
sector must satisfy an analogous infinite set of sum-rules. In this case, the sum-rules can be 
solved, and there is a large body of work successfully comparing the solutions to numerical 



results from overlap and related fermions (see, for example, the review in Ref. (22|). It would 
be of considerable interest if the solutions could be extended to the sum-rules discussed here. 

One might wonder what can be learned about the other scenario predicted by xPT — 
that involving a first-order transition. The answer appears to be very little. In this case 
p(A) always has a gap, there is no SSB, there are no massless Goldstone pions, and thus no 
e-regime. One expects from xPT or from the p.d.f. analysis that (Cq) = (C^) = for both 
a = 0, 3. The consistency of these results with the expressions in terms of eigenvalues is 
almost trivial, because with a gap there are no infrared divergences. 

One might also wonder what happens to the present analysis in the continuum limit. In 
particular, how does it connect with that of Leutwyler and Smilga? The short answer is that 



I thank Vicente Azcoiti and collaborators for stressing to me the importance of the presence of this infinite 
tower of sum-rules. 
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there is no direct connection, since the hmits a ^ and A^sitc — > oo do not commute when 
in the Aoki-phase. If one takes iVsite — oo first, as I have throughout, one is always in the 
p-regime for the neutral pion, so that the vacuum manifold is SU{2)/U{1), while, with a ^ 
first (and at the same time sending m — > so as to remain in the Aoki-phase), the manifold 
is SU{2). Another way of seeing the difference is to note that, if NsHc — ^ before a — 0, 
then would-be zero modes are completely buried in the continuum of near- zero modes, while 
if a — > first then the would-be zero modes lie below the continuum of near-zero modes. 

One can, nevertheless, ask what happens if a — first. One still expects (Cq)^ = (from 
XPT or from the properties of p.d.f.s), and must understand this result. The answer turns 
out to be that zero-mode contributions to both (CQ)disc,e and (CQ)conn,e conspire to cancel the 
contribution from the continuum of near-zero modes to (Cg)conn,e- The cancellation occurs 
as long as one of the sum-rules of Leutwyler and Smilga holds. Thus one could obtain this 
sum-rule by enforcing (Cq)^ = 0. In fact, it is possible to obtain the whole set of higher-order 
sum-rules by enforcing the xPT relations between condensates of higher powers. This shows 
the close relation between the methods used here and those of Ref . 1^ . 
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APPENDIX A: SPECTRAL ASYMMETRY AND TOPOLOGICAL SUSCEPTI- 
BILITY AT a / 

This appendix provides a more detailed argument for the result (15. 6p used in the main 
text for the spectral asymmetry of H^/. The argument consists of two parts. The first aims 
to justify the approximate proportionahty 

(mo)n oc (Qtop("^o)') = Vxtimo) , (Al) 

where Xt is the topological susceptibility, and the arguments mo indicate that all quantities 
are evaluated in the Aoki-phase. This proportionality is for fixed mo, and thus concerns 
the dependence on a and V. The second part of the argument gives a derivation of the 
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parametric form of the Xt in the Aoki-phase using "Wilson fermion xPT" (W^PT) j4j. |23|. 
Combining these two parts leads to the desired result. 

From the work of Refs. 16, 1^, 18 1 we have a fairly clear picture of how the spectral 



asymmetry occurs. Given a configuration, one considers the spectrum of low-lying eigenval- 
ues of the "valence" Hermitian Wilson-Dirac operator, Hw{my), in which mg is replaced by 
the bare valence quark-mass my. One studies this spectrum as a function of my, starting 
at positive values and decreasing to my- pa —1. Based on the numerical results of Ref. |l8| . 



and subsequent theoretical work 20|, |2J] , I assume that there is a "valence Aoki-phase" for 



a region of negative my. For positive my the spectral asymmetry vanishes, but, once my is 
negative and one enters the "supercritical region", eigenvalues can cross zero. This means 



that, by the time one reaches the Aoki-phase, some spectra! 



asymmetry can have built up, 



and numerical results indicate that this in fact happens [l8j. The asymmetry increases as 
one moves through the Aoki-phase, and becomes almost independent of my shortly after 
leaving the Aoki-phase. This approximate independence continues down to my ^ —1, and 
is due to a cancellation of modes crossing in both directions rather than an absence of cross- 
ings. The resulting spectral asymmetry for my ^ —1 is known to provide a robust definition 



of the topological charge o 
From the results of Ref. 



the configuration 



id. 



18| , this behavior appears to hold for a wide variety of ensembles. 



both quenched and unquenched. I assume here that it also holds for unquenched ensembles 
in which mo is in range leading to a dynamical Aoki-phase. If so, then, for a given mo, ?T-asym 
will approximately track Qtop as a and V are varied, leading to the result ( lAip . Note that 
the proportionality constant between these two quantities will depend on mo, i.e. on where 
one lies in the Aoki-phase. Indeed, the behavior described in the previous paragraph implies 
that nasym/Qtop is small near the upper boundary of the Aoki-phase (larger mo) and close 
to unity at the lower boundary. 

Further justification for the assumed proportionality comes from considering the physical 
extent of the eigenmodes at zero-crossing. According to Ref. [l^, this typically decreases 
as my is decreased. Crudely speaking, one can think of the crossings which occur before 
entering and within the Aoki-phase as corresponding to "lumps" of topological charge that 
extend over many lattice spacings, and survive the continuum limit. By contrast, those 
that occur after traversing the Aoki-phase are mostly lumps of size ~ a. Thus, if my 
is within the Aoki-phase, the spectral asymmetry gives an approximate measure of that 
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part of Qtop resulting from lumps of greater than some minimal size. Since the topological 



susceptibility in the continuum limit (appropriately regularized 25|) is determined by lumps 
of size ~ 1/Aqcd, it is plausible that this "truncated" or "coarse-grained" topological charge 
should lead to a susceptibility proportional to the exact result, and thus to Eq. flAip . 

1 now move to the second part of the argument. Accepting (lAip . the next task is to 
determine the expected dependence of Xt in the Aoki-phase on a (and possibly V). Given 
some assumptions, this can be done using W^PT, generalizing the standard continuum 



xPT analysis |10l. |26[|. As a byproduct, the form of the dependence of Xt on mo will also be 
obtained, but this does not carry over to (|?T-asym('^v)P) because the proportionality constant 
in (lAip depends on ttiq. 

In the continuum analysis, one introduces the OFF term into the action, rotates it into 
the quark mass matrix using an anomalous singlet axial rotation, and then evaluates Z{9) 
using xPT. At leading order, and in the p-regime at large volume, Z{9) = exp[— V^V^jjj(6')], 
where Vj^j^ is the (6'-dependent) minimum of the potential in the chiral Lagrangian. Then 
one has 



{dey 



(A2) 



y-oo V {dey 

To generalize this to include discretization effects, one must certainly include O(a^) correc- 
tions to the potential, for these contribute to the vacuum energy at leading order in W^PT. 
It is less clear, however, how to introduce 6 into W^PT. One might consider adding a bare 
6FF term to the lattice action, mapping this to the Symanzik action, and then into the 
chiral Lagrangian — i.e. following the standard steps in W^PT. Such a bare operator will, 
however, mix with lower [ijj'-f^ilj) and higher {ip'j50'^i,F^utlj) dimension operators, leading re- 



spectively to 0{l/a) effects that need to be subtracted non-perturbatively 27|] and 0{a) 
discretization errors. 

I think, however, that these complications do not occur here because the asymmetry 
that appears in Eq. (1A1|) is proportional to a Qtop that is regulated in the ultraviolet. 
This precludes mixing with lower-dimensional operators. In other words, the definition 
Qtop = '^asym("^v' ~ — 1) gives a result that requires no subtractions. Qtop will, however, have 
discretization errors. For example, using an improved Wilson-Dirac operator in the valence 
H\Y would lead to an different assignment of topological charge on some configurations. 
What I assume here is that this is an O(a^) error rather than an 0(a) one, because one is. 
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in effect, using valence overlap fermions to measure Qtop, and overlap fermions have only 
O(a^) errors. 

What this discussion leads to is the assumption that, in order to calculate Xt, the appro- 
priate potential to use in the Symanzik continuum effective action is 

Vsymanzik ~ ^(^ + i^-i^n)^) + a#cT ■ + {ij^jf + lOFF + 0(^0^) . (A3) 

Here I use a schematic notation in which all constants of 0(1) have been dropped, and 
I have only shown one of the possible forms of the chiral-symmetry-breaking four-fermion 
operators. I have also dropped terms of higher order in the chiral counting m ~ a^. In 
words, my assumption is that the lattice Qtop matches onto / FF + 0{a?) in flA3l) . The 
0{9a?) terms will give rise to mass independent 0{o?) corrections to Xt-, which I will drop 
for now but restore at the end. 

Given Eq. ( 1A3I) [minus the 0{9a?') terms] the remainder of the analysis is straightforward. 
One first rotates 9 into the fermionic terms, with each (scalar or pseudoscalar) bilinear pick- 
ing up a factor of exp(z6'75/2).^^ One then matches to the chiral effective theory, obtaining 
the same form as for ^ = except for the substitution E e^*^S, where 9 = 9/2. Thus the 
potential becomes 



+ ^ |2Tr 
16 1 



e-'^E + e'^E'^y - Tr(e-^^S + e*^S^ 



Tr(e"'^^S + e'^S^) 

2~1 



2 



(A4) 



Here I have absorbed the 0{a) term into the mass matrix M = m + ifir^ in the usual way, 
which remains possible even when 9^0. Compared to Eq. (12.11) in the main text, there 
is an additional term, that proportional to ~ a^Aq^j-,. This term vanishes when 9 = 
due to the properties of SU (2) matrices, but is non-vanishing when 9^0. Recall also that 
2/2^0 = ci/m. 

To determine Xt using (lA2p one must minimize with respect to variations in E. This 
assumes fi is kept non-zero while F ^ oo so that we are in the p-regime. Inserting S = 
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It is possible that one can avoid the extended discussion given above by starting with 9 inserted into the 



quark mass matrix of the lattice theory, following the work of Ref. 28|. I have not, however, pursued to 
resolution the issues of renormalization that arise in this approach. 
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exp(i^o^ ■ ^) one finds, up to an irrelevant constant, that 

V = — ci cos 6^0 cos 6' + ci{fi/m) sin 6*07^3 cos 9 + C2 cos^ 6*0 cos^ 9 + 02(1 — cos^ ^o)(l — cos^ 9) . 

(A5) 

Minimizing with respect to 9q, and considering only infinitesimal /i, leads to 

V„,.(9) = + «-fHl:ii£lZ?f!fl + 0{g') , (A7) 

4c2 4 

and thus 

X. = ''^^'-^l'/''^^'^ + 0(a^A«QCD) = + O(a^A^QCD) • (A8) 

One finds that the "calculable part" of Xt is determined by the new low-energy constant, 
I have also reinserted the O(a^) term resulting from the 0{9a^) contributions to Eq. (]A3p . 
The O(a^) term is of the same size as the contribution, but does not depend on mo, being 
simply a discretization error in Qtop and not related to the alignment of the vacuum. The 
overall conclusion is thus that Xt ~ ct^^QCD) so that on a typical configuration, |nasym| ~ 
y/Vxt ~ Aq(^£,. This is the result used in the main text. 

Finally, I address what happens in the e-regime for the charged pions. Here one must 
include in the calculation of Z{9) the integral over the direction, n, of the condensate in 



the SU{2)/U{1) manifold. Following the method of Ref. lOj I find that this adds to Vxt 
a contribution proportional to [fiV Aqq^))"^ . In order to be in the e-regime, however, this 
contribution must have magnitude much smaller than unity. Thus it has no impact on the 
topological susceptibility when V 00. 



APPENDIX B: ALTERNATIVE FORMULATION OF CONSISTENCY CONDI- 
TIONS 

In this Appendix I describe an alternative, and arguably more natural, formulation of 
the conditions which must be satisfied by eigenvalue distributions in order that the standard 
Aoki-phase scenario remain valid. The formulation uses the microscopic spectral density, and 
the corresponding higher-order eigenvalue correlations. For the continuum Dirac operator, 
these are the quantities whose properties are universal in QCD-like theories, and governed 
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by random matrix theory. To define them one "zooms in" on the region of eigenvalues of 



size A < (aAQCD)/A^. 



site' 



.12 



Pi(^)=P(A = ^). (Bl) 

* si 



'site 



Thus p\{x)dx is the average total number of eigenvalues between A = x/N^ite and {x + 
dx)/Nsitf.. When expressed in terms of x, eigenvalues on a single configuration form a 
discrete set of levels, spaced by ~ l/p(0). One obtains a continuous distribution only 
after averaging over configurations, and the resulting distribution has structure (oscillations 
about p(0)) which is the remnant of the discrete levels. This is different from p(A), which, as 
noted in the main text, becomes a continuous function pc/(A) on a single configuration when 
A'^site —>■ oo- Furthermore, p(A) does not display the oscillations for small A seen in pl{x), 
since they get averaged out when NsHf. oo. The essential point is that pl{x) contains 
extra information about the IR region that is lost in p(A). 

I assume in the following that the lattice quark mass is chosen so that p(0) 7^ 0, implying 
that one is in the Aoki-phase in the standard xPT description. The eigenvalues of interest, 
which I call the "IR eigenvalues", are those for which p(A) is approximately constant, i.e. for 
which higher-order chiral corrections are small. This requires |A| <^ oAqcd, which translates 
into a maximum x of magnitude Xmax = c(aAQCD)A'site, with c ^ 1 a positive constant. Since 
pI{x) is on average a constant, the number of eigenvalues in the range — x^ax < a; < x^ax is 
/-^max dxpl{x) ~ 2a;maxP(0), and in particular is proportional to Agite- 

The distribution of the IR eigenvalues is encoded by pl{x), along with higher order corre- 
lation functions, pl(x-\, . . . ,Xk)- The latter are standard quantities and I use the definitions 



given in Ref. 29|, that, in particular, do not include the subtraction of the "connected 
part" . Thus, if the eigenvalues were completely uncorrelated, one would have, for example, 
p|(xi,X2) = pI{xi)pI{x2) — S{xi — X2)pI{xi). The second term is present because p| is a cor- 
relation between the eigenvalues of distinct eigenvectors. To simplify some of the following 
formulae, I also use correlation functions in which these kinematic correlations are removed, 
e.g. 

p2(Xi,X2) = p2{Xi,X2) + S{Xi - X2)pliXi) . (B2) 



Other scaling factors are also used in the continuum literature, e.g. including a factor of the condensate, 
S oc p(0), so that the eigenvalue spacing is of 0(1). I prefer not to do this since p(0) is not a constant, 
but rather depends on the position in the Aoki-phase. Note also that some authors define p{X) without 
dividing by A'^gite, in which case an additional factor ofl/NsHe is needed on the right-hand-side of eq. (|Bip . 
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Then for uncorrelated eigenvalues one has . . . , Xk) = ni=i,fc Pii^i) for all k. Both the 

pI and pI are symmetric under interchange of any two arguments. It it also useful to note 
the volume dependence of the normalization of the pf., 

n dxi\ pk{xi,..., Xk) oc (A^site) , (B3) 



Ki=l,k 



which is consistent with the expectation that p| approaches a constant when all arguments 
have magnitudes much larger than l/p(0). 

In the continuum, the corresponding correlators are even functions of each Xi separately, 
allowing one to work only with Xi > 0. On the lattice, however, parity invariance of Hw 
implies only that the pf, do not change when all arguments change sign simultaneously, e.g. 

plixi, X2, . . . , Xk) = pli-xi, -X2, . . . , -Xk) ■ (B4) 

The same holds for the p^. Thus, antisymmetric parts such as 

Ap'^ixi, X2) = [P2ixi, X2) - P2{~Xi,X2)] /2 = -Ap^(-a;i, xs) (B5) 

need not vanish, unlike in the continuum. In fact, this particular quantity must be non- 
vanishing since it encodes the spectral asymmetry, which itself is non-zero (as discussed in 
the main text and in Appendix : 

(nlyj = {iN+ - Ar.)2) = A dx^dx2 Ap',{x,, X2) , (B6) 

(Note that the integrals here are over positive Xi only.) We do, however, learn from Eq. (]B6P 
that Ap2 cannot tend to a constant for large \xi\, unlike P2 itself. This is because the inte- 
gration area grows as N^^^^ while the integral grows only as (?T-asym) ^site [from Eq. (1A1I) ]. 
One possible behavior is that Ap^ falls off for large Xi, and this indeed is what is suggested 
by the discussion of the spectral asymmetry in the main text. 

It is straightforward to convert the sum-rules given in the main text into constraints on 
the pI and/or pi. Results are simplified by defining (following continuum usage) 

So = |^sin^o/'So = 7rp(0). (B7) 

Zjp 



The sum-rule ( 15.21) then becomes 

max ^ 
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This has the same form as the first Leutwyler-Smilga sum-rule for QCD [10|, except that 
2/3 is replaced in QCD by 1/2(2+ |z/|), with v the topological charge. One must also exclude 
exact zero-modes from p\ in the continuum. Since p\[x) is symmetric, one could restrict 
the integral to positive values only (and divide the right-hand-side by two) — this is how the 
sum-rules are usually expressed in the continuum. 

In this and subsequent sum-rules an implicit limit of N^xtc oo has been taken. 
This removes the UV contribution to the integral (arising from the large x behavior, 
pI oc (x/iVgite)'^), so that the dominant contribution to the integral is from the IR region. 
Note also that the precise upper limit, Xmax, is then irrelevant. 

The next sum- rule is obtained from Eq. (15.51) . and is 

dx.dx/-^^^ = 4 dx.dx,^^^^ = f . (B9) 

J — Xma^L X1X2 Jo X1X2 3 

The oddness of the integrand in both xi and X2 picks out the antisymmetric part Apg. This 
means that this sum-rule has no continuum analog. 

In the discussion given in Sec. IVj the consistency of the Aoki phase required what 
might appear to be an artificial construct, namely an antisymmetry in pu{^), subleading as 
iVsite — > 00, and yet contributing to the sum- rules due to its IR divergence. By contrast, the 
consistency conditions seem quite natural in the present formulation. The lattice symmetries 
allow a new quantity to be present, i.e. A/92(xi, X2), a function which one expects to remain 
non- vanishing in the NsHc 00 limit. Just as pi needs to satisfy a consistency condition 
[eq. IB8j . it seems natural that the new function should too. Certainly from a mathematical 
point of view there should be no barrier to satisfying (IBSp . since one has a function to play 
with and the only other constraint is eq. flB6|) . One solution is the second example given in 
Sec. lYl in which Ap2 oc l/^/XlX2^^ 

The pattern of one sum-rule for each new function continues at higher order. The sum 
rules are obtained from enforcing 

(O. = (C^rCo"')e = and (C^r)e = ^^((C^3)p)'" (BIO) 

for integer values of n, n', results that can be obtained from xPT or using probability 
distribution functions. In particular, combining the sum-rules for Cg, CqC^ and C3, one 



This form can hold only for large xi and X2. In particular, oddness in a;i and X2 separately implies that 
Ap| vanishes when xi — > or a;2 ^ 0, properties also needed if the integral is to converge in the IR. 
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finds 



4 /-n.ax PI{XI,X2) 



_Et = / dx,dx2'-^^-^ . (Bll) 



— S'^ = / dxidx2dx^ ^^^^^'' '^'^''^ ^ 

15 J-XmsL^ 



= / cixidx2c/x3 ^^^:^:,-^!r^^ ■ (B12) 



-E^ = f dxidx2dx^dxi ^'^^ — — — — — — . (B13) 

Note that it is simpler here to use the than the p|. The first constraint picks out 
the continuum-hke, symmetric part of p|, and indeed a result of similar form holds in the 
continuum. The other two constraints pick out parts of the three- and four-point correlators 
that vanish in the continuum, but need not vanish on the lattice. These are new, essentially 
independent, functions, and there is no barrier that I can see to their satisfying the new 
sum-rules. 

In principle, one can continue this procedure to arbitrarily high order, and obtain an 
infinite set of sum-rules. This set is analogous to those that hold for each value of v in the 
continuum, except that on the lattice one has the additional functions (those odd under sign 
fiips of a subset of the Xi) and corresponding additional sum-rules. In the continuum, these 
sum-rules have been solved to obtain the p|., and from them the distribution of individual 
eigenvalues (as reviewed in Ref. 221]). It is an interesting challenge to extend this analysis 
to the lattice theory. 
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